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Abstract 



The paper presents the short review of our to-day knowledge of vacuum con- 
densates in QCD. The condensates are defined as vacuum averages of the operators 
which arise due to nonperturbative effects. The important role of condensates in de- 
termination of physical properties of hadrons and of their low-energy interactions in 
QCD is underlined. The special value of quark condensate, connected with the exis- 
tence of baryon masses is mentioned. Vacuum condensates induced by external fields 
are discussed. QCD at low energy is checked on the basis of the data on hadronic 
T-decay. In the theoretical analysis the terms of perturbation theory (PT) up to 
are accounted, in the operator product expansion (OPE) - those up to dimension 
8. The total probability of the decay r — > hadrons (with zero strangeness) and of 
the r-decay structure functions are best described at a s {m 2 ) = 0.330 ± 0.025. It is 
shown that the Borel sum rules for r-decay structure functions along the rays in the 
g 2 -complex plane are in agreement with the experiment with the accuracy ~ 2% at 
the values of the Borel parameter \M 2 \ > 0.8GeV 2 . The magnitudes of dimensions 
6 and 8 condensates were found and the limitations on gluonic condensates were ob- 
tained. The sum rules for the charmed quark vector currents polarization operator 
was analysed in 3 loops (i.e., in order a 2 ). The value of charmed quark mass (in 
MS regularization scheme) was found to be: m c (m c ) = 1.275 ± O.OlSGeV^ and the 
value of gluonic condensate was estimated: (0|(a s /7r)C7 2 |0) = 0.009 ± 0.007GeV 4 . 
The general conclusion is: QCD described by PT + OPE is in a good agreement 
with experiment at Q 2 > lGeV 2 . 
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1. A few words about Igor' Vasilyevich. 



It is a great honour and at the same time a great pleasure for me to write a 
paper to the issue of Physics of Atomic Nuclei dedicated to 100-anniversary of Igor' 
Vasilyevich Kurchatov birthday. Kurchatov was a very extraordinary person: an 
organizer of the highest class, I know nobody with such excellent organization abil- 
ities. Without him, the soviet atomic project would not be, perhaps, realized, at 
the least in such short time. Igor' Vasilyevich had the strongest sense of respon- 
sibility not only for the work entrusted to him - the atomic project - but for the 
much wider - for the fortune of the science in our country and moreover, for the 
fortune of the whole mankind. I shall give an episod not well known. As witnesses 
A.P.Aleksandrov [1], Kurchatov was deeply depressed when coming back from the 
tests of the first hydrogen bomb (those who were present at the tests noticed the 
same). He said: "What a terrible thing we have made. The only item we should 
bother about, is to forbid all of this and to exclude nuclear war". In March 1954 
Kurchatov, Alikhanov, Kikoin and Vinogradov had written the paper, where they 
concluded, that "- the mankind is in front of the menace of the end of all of the life 
on the Earth". The paper was also signed by the Minister of the Middle Machine 
Building V.A.Malyshev who had sent it to Malenkov, Khrushchev and Molotov. 
Khrushchev, however, had rejected the paper, calling the words on a possible ruin 
of the world civilization "theoretically wrong and politically harmful" (see [2]). The 
position of the soviet leaders remained as before: the world war should lead to the 
ruin of the capitalism. 

The same responsibility was inherent to Igor' Vasilyevich when constructing 
atomic reactors and atomic powerstations - I gave examples of this earlier [3]. I 
think that if Kurchatov would alive, RBMK reactors, principally not safe as physical 
systems, would not be constructed and we would avoid the Chernobyl catastrophe. 
But on the other hand, Igor' Vasilyevich was a person of his time ... (see [3]). 

He liked the science and, in the first turn, his main speciality-nuclear physics. He 
was deeply interested in the development of the elementary particle physics and he 
thought that it is necessary to develop such investigations in the USSR. He supported 
the suggestion of Alikhanov and Vladimirsky of construction at ITEP of the 7 GeV 
hard-focusing proton accelerator, and then, using the ITEP project, of the 50 GeV 
proton accelerator (later, 70 GeV) near Serpukhov. In 1954 such decision had been 
adopted at the meeting of the chaired by Kurchatov Scientific- Technical Council of 
the Middle Machine-Building Ministry. 

2. Introduction 

Nowadays, it is reliably established that the true (microscopic) theory of strong 
interaction is quantum chromodynamics (QCD), the gauge theory of interacting 
quarks and gluons. It is also established, that unlike, e.g., quantum electrodynam- 
ics (QED), the vacuum in QCD has a nontrivial structure: due to nonperturbative 
effects, i.e. not admitting the expansion in the interaction constant (even if it is 
small) in QCD vacuum persist non-zero fluctuations of gluonic and quark fields. 
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(Examples of such kind of nonperturbative fields are instantons [4] - classical so- 
lutions of equations for gluonic field, which realize the minimum of actions in the 
QCD LagrangianQ The nontrivial structure in QCD manifests itself in the pres- 
ence of vacuum condensates, analogous to those in the condensed matter physics 
(for instance, spontaneous magnetization). Vacuum condensates are very impor- 
tant in elucidation of the QCD structure and in description of hadron properties at 
low energies. Condensates, in particular, quark and gluonic ones, were investigated 
starting from 70-ties. Here, first, it should be noted the QCD sum rule method by 
Shifman, Vainshtein, and Zakharov [5], which was based on the idea of the leading 
role of condensates in the calculation of masses of the low-lying hadronic states. In 
the papers of 70-80-ies it was adopted that the perturbaive interaction constant is 
comparatively small (e.g., a s {lGeV) ~ 0.3), so that it is enough to restrict oneself 
by the first-order terms in a s and sometimes even disregard perturbative effects in 
the region of masses larger than 1 GeV. At present it is clear that a s is consider- 
ably larger (a s (lGeV) ~ 0.6). In a number of cases there appeared the results of 
perturbative calculations in order a 2 s and a^. New, more precise experimental data 
at low energies had been obtained. Thereby, on one hand, it is necessary, and on 
the other it appears to be possible to compare QCD with the experiment in the low 
energy region on a higher level of precision. The results of such a comparison are 
presented in this paper. 

In Section 3 I define condensates, describe their properties and give numerical 
values which were obtained previously. In Section 4, the data on hadronic decays 
of r-lepton are compared with theoretical expectations obtained on the basis of the 
operator product expansion in QCD with the account of perturbative terms up to 
a^. The values of condensates and the coupling constant a s (m 2 ) are obtained. In 
Section 5, polarization operator of the vector current of charmed quarks is analysed 
in three-loop approximation (i.e., with the account of the terms ~ a 2 ), the value of 
the charmed quark and the value of gluonic condensate are found. Section 6 presents 
our conclusions. 

3. Definition of condensates, their main properties 

In QCD (or in a more general case, in quantum field theory) by condensates there 
are called the vacuum mean values (0|Oj|0) of the local (i.e. taken at a single point of 
space-time) of the operators Oj(x), which arise due to nonperturbative effects. The 
latter point is very important and needs clarification. When determining vacuum 
condensates one implies the averaging only over nonperturbative fluctuations. If for 
some operator Oi the non-zero vacuum mean value appears also in the perturbation 
theory, it should not be taken into account in determination of the condensate - in 
other words, when determining condensates the perturbative vacuum mean values 
should be subtracted in calculation of the vacuum averages. One more specifica- 
tion is necessary. The perturbation theory series in QCD are asymptotic series. So, 

1 This statement refers to Euclidean space, in the Minkowsky space, instantons realize tunneling 
transitions between Hilbcrt spaces with different topological quanum numbers. 
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vacuum mean operator values may appear due to one or another summing of asymp- 
totic series. The vacuum mean values of such kind are commonly to be referred to 
vacuum condensates. 

Separation of perturbative and nonperturbative contributon into vacuum mean 
values has some arbitrariness. Usually [6,7], this arbitrariness is avoided by intro- 
duction of some normalization point /i 2 (/i 2 ~ lGeV 2 ). Integration over momenta 
of virtual quarks and gluons in the region below fi 2 is referred to condensates, above 
/i 2 - to perturbative theory. In such a formulation condensates depend on the nor- 
malization point ji\ (0|Oj|0) = (0|Oj|0) M . Other methods for determination of 
condensates are also possible (see below). 

In perturbation theory, there appear corrections to the condensates as a series 
in the coupling constant a s (/i): 

oo 

(0|O i |0>g= (0|O,|0> M ]T C®{Q,fia n M (1) 

n=0 

The running coupling constant a s at the right-hand part of (1) is normalized at the 
point fi. The left-hand part of (1) represents the value of the condensate normalized 
at the point Q. Coefficients C^(Q,fj) may have logarithms lnQ 2 /[i 2 in powers up 
to n for . Summing up of the terms with highest powers of logarithms leads to 
appearance of the so-called anomalous dimension of operators, so that in general 
form it can be written 

(oio.io) = (°I^I°)m(^) 7 £ <$>(Q,nWM, (2) 

where 7 - is anomalous dimension (number), and has already no leading log- 
arithms. If there exist several operators of the given (canonical) dimension, then 
their mixing is possible in perturbation theory. Then the relations (1),(2) become 
matrix ones. 

In their physical properties condensates in QCD have much in common with 
condensates appearing in condensed matter physics: such as superfluid liquid (Bose- 
condensate) in liquid A He, Cooper pair condensate in superconductor, spontaneous 
magnetizaion in magnetic etc. That is why, analogously to effects in the physics of 
condensed matter, it can be expected that if one considers QCD at finite temperature 
T, with T increasing at some T = T C there will be phase transition and condensates 
(or a part of them) will be destroyed. Particularly, such a phenomenon must hold 
for condensates responsible for spontaneous symmetry breaking - at T = T c they 
should vanish and symmetry must be restored. (In principle, surely, QCD may have 
a few phase transition). 

Condensates in QCD are divided into two types: conserving and violating chi- 
rality. As is known, the masses of light quarks u, d, s in the QCD Lagrangian are 
small comparing with the characteristic scale of hadronic masses M ~ lGeV. In 
neglecting light quark masses the QCD Lagrangian becomes chiral-invariant: left- 
hand and right-hand (in chirality) light quarks do not interact with each other, 
both vector and axial currents are conserved (except for flavour-singlet axial cur- 
rent, non-conservation of which is due to anomaly). The accuracy of light quark 
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masses neglect corresponds to the accuracy of isotopical symmetry, i.e. a few per 
cent in the case of u and d quarks and of the accuracy of SU(3) symmetry, i.e. 
10-15 % in the case of s-quarks. In the case of condensates violating chiral symme- 
try, perturbative vacuum mean values are proportional to light quark masses and 
are zero within m u = rrid = m s = 0. So, such condensates are determined in the 
theory much better than those conserving chirality and, in principle, may be found 
experimentally with higher accuracy. 

Among chiral symmetry violating condensates of the most importance is the 
quark condensate (0|gg|0) (q = u,d are the fields of u and d quarks). (0|gg|0) may 
be written in the form 

(0\qq\0)^(0\q L q R + q R q L \0) (3) 

where q^, q R are the fields of left-hand and right-hand (in chirality) quarks. As 
follows from (3), the non-zero value of quark condensate means the transition of 
left-hand quark fields into right-hand ones and its not a small value would mean to 
chiral symmetry violation in QCD. (If chiral symmetry is not violated, then at small 
m u ,md (0|gg|0) ~ m u ,md). By virtue of isotopical invariance 

(0|u«|0) = (0\dd\0) (4) 

For quark condensate there holds the Gell-Mann-Oakes-Renner relation [8] 

1 m 2 f 2 
z m u + ma 

Here m n ,f n are the mass and constant of 7r + -meson decay (m„ = 140MeV, f n = 
131MeV"), m u and are the masses of u and d-quarks. Relation (5) is obtained 
in the first order in m u ,md,m s (for its derivation see, e.g. [9]). To estimate the 
value of quark condensate one may use the values of quark masses m u = 4.2MeV, 
m d = 7.5MeV [9]. (These values were suggested by Weinberg [10], within the errors 
they coincide with other estimates - see, for example, [11]). Substituting these 
values into (5) we get 

(0\qq\0) = -(243Me\/) 3 (6) 

The value (6) has characteristic hadronic scale. This shows that chiral symmetry 
which is fulfilled with a good accuracy in the light quark lagrangian (m u ,m^/M ~ 
0.01, M - is hadronic mass scale, M ~ 0.5 — lGeV), is spontaneously violated on 
hadronic state spectrum. 

An other argument in the favour of spontaneous violation of chiral symmetry 
in QCD is the existence of massive baryons. Indeed, in the chiral-symmetrical the- 
ory all fermionic states should be either massless or parity-degenerated. Obviously, 
baryons, in particular, nucleon do not possess this property. It can be shown [12, 9], 
that both these phenomena - the presence of the chiral symmetry violating quark 
condensate and the existence of massive baryons are closely connected with each 
other. According to the Goldstone theorem, the spontaneous symmetry violation 
leads to appearance of massless particles in the physical state spectrum - of Gold- 
stone bosons. In QCD Goldstone bosons can be identified with a 7r-meson triplet 
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within m U) md — ► 0, m s ^ (SU(2)-symmetry) or with an octet of pseudoscalar 
mesons (ir, K, rj) within the limit m u , m^, m s — > (SU(3)-symmetry). The presence 
of Goldstone bosons in QCD makes it possible to formulate the low-energy chiral 
effective theory of strong interactions (see reviews [13], [14], [9]). 

Quark condensate may be considered as an order parameter in QCD correspond- 
ing to spontaneous violation of the chiral symmetry. At the temperature of restora- 
tion of the chiral symmetry T = T c it must vanish. The investigation of the tem- 
perature dependence of quark condensate in chiral effective theory [15] (see also 
the review [9]) shows that (0|gg|0) vanishes at T = T c 150 — 200MeV. Similar 
indications were obtained also in the lattice calculations [16]. 

Thus, the quark condensate: l)has the lowest dimensions (d=3) as compared 
with other condensates in QCD; 2) determines masses of usual (nonstrange) baryons; 
3) is the order parameter in the phase transition between the phases of violated and 
restored chiral symmetry. These three facts determine its important role in the 
low-energy hadronic physics. 

Let us estimate the accuracy of numerical value of (6). The Gell-Mann-Oakes- 
Renner relation is derived up to correction terms linear in quark masses. In the 
chiral effective theory one succeeds in estimating the correction terms and, thereby, 
the accuracy of equation (5) appears of order 10%. However, it is not a single origin 
of errors in determination of quark condensate value. The quark condensate, as well 
as quark masses depend on the normalization point and have anomalous dimensions 
equalling to 7 m = ~lq q — f ■ I n the mass values taken above the normalization 
point /i was not fixed exactly (in fact, it was taken // ~ lGeV). In addition, the 
accuracy of the above taken value m u + m d = 11.7 MeV which enters (5) seems to 
be of order 10 — 20%. The value of the quark condensate may be also found from 
the sum rules for proton mass. The analysis made [17] gave for it a value very close 
to (6) (with the 3% difference) at the normalization point fi = lGeV. The accuracy 
of these sum rules seems to be of order 10-15%. Concludingly, it may be believed, 
that the value of the quark condensate is given by (6) at the normalization point 
\i = lGeV with the 10-20% accuracy. The quark condensate of strange quarks is 
somewhat different from (0|ww|0). In [12] it was obtained 

(0|ss|0)/(0|uu|0) = 0.8 ±0.1 (7) 

The next in dimension (d = 5) condensate which violates chiral symmetry is quark 
gluonic one: 

- g(0\qa, uY G; u q\0) =ml(0\qq\) (8) 

Here G™ u - is the gluonic field strength tensor, A™ - are the Gell-Mann matrices, 
®p.v = (if (l/ilv ~ Ivln- The value of the parameter wig was found in [18] from the 
sum rules for baryonic resonances 

ml = 0.8GeV 2 (9) 

Consider now condensates conserving chirality. Of fundamental role here is the 
gluonic condensate of the lowest dimension: 

<0|-r<MJ0> (10) 
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Due to that the gluonic condensate is proportional to the vacuum mean value of 
the trace of the energy-momentum tensor 9^ u its anomalous dimension is zero. The 
existence of gluonic condensate had been first indicated by Shifman, Vainshtein, 
and Zakharov [5]. They had also obtained from the sum rules for charmonium its 
numerical value: 

ry 

(0\fG; u G; u \0) = 0M2GeV 4 (11) 

As was shown by the same authors, the nonzero and positive value of gluonic con- 
densate mean, that the vacuum energy is negative in QCD: vacuum energy density 
in QCD is given by e = — (9/32)<0| (a s /7r)G 2 |0). Therefore, if quark is embedded 
into vacuum, this results in its excitation, i.e, in increasing of energy. In this way, 
the explanation of the bag model could be obtained in QCD: in the domain around 
quark there appears an excess of energy, which is treated as the energy density B 
in the bag model. (Although, the magnitude of B, does not, probably, agree with 
the value of e which follows from (11)). In ref.[5] perturbative effects were taken 
into account only in the order a s , the value for a s being taken twice as smaller as 
the modern one. Later many attempts were made to determine the value of gluonic 
condensate by studying various processes and by applying various methods. But the 
results of different approaches were inconsistent with each other and with (11) and 
sometimes the difference was even very large - the values of condensate appeared to 
be by a few times larger. All of this needs reanalysation of (0| — G 2 |0) determination 
basing on contemporary values which will be done in Sections 4,5. 
The d=6 gluonic condensate is of the form 

g ^f^(0\Gl u Gl x G^\0), (12) 

(f abc - are structure constants of SU(3) group). There are no reliable methods to 
determine it from experimental data. There is only an estimate which follows from 
the method of deluted instanton gas [19]: 

9 3 f abc (0\G%G b v\G^\0) = ^(12^)1(01^10), (13) 

where p c is the instanton effective radius in the given model (for estimaion one may 
take p c ~ (1/3- 1/2) fm). 

The general form of d=6 condensates is as follows: 

a s (0\qiO a qi ■ q k O a q k \0) (14) 

where qi,q k are quark fields of u,d,s quarks, O a - are Dirac and SU(3) matrices. 
Following [5], Eq.(14) is usually factorized: in the sum over intermediate state in 
all channels (i.e, if necessary, after Fierz-transformation) only vacuum state is taken 
into account. The accuracy of such approximation ~ l/A" 2 , where N c is the number 
of colours i.e.~ 10%. After factorizaion Eq.(14) reduces to 

a s (0\qq\0) 2 (15) 
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The anomalous dimension of (15) is 1/9 and it can be approximately put to be zero. 
And finally, d=8 quark condensates assuming factorization reduce to 

a a (0\qq\0)-m*(0\qq\0) (16) 

(The notaion of (8) is used). It should be noted, however, that the factorization 
procedure in the d=8 condensate case is uncertain. For this reason, it is necessary 
to require their contribution to be small. 

Let us also dwell on one more type of condensates - those, induced by external 
fields. The meaning of such condensates can be easily understood by comparing with 
analogous phenomena in the physics of condenced media. If the above considered 
condensates can be compared, for instance with ferromagnetics, where magnetiza- 
tion is present even in the absence of external magnetic field, condensates induced 
by external field are similar to dia- or paramagnetics. Consider the case of the con- 
stant external electromagnetic field F^ v . In its presence there appears a condensate 
induced by external field (in the linear approximation in F^ v ): 

(0|g<vg|0) F = e ?X ^(0|^|0) (17) 

As was shown in ref. [20], in a good approximation (0|go" M ^g|0)ir is proportional to e q 
- the charge of quark q. Induced, by the field vacuum expectation value (Olga^glO) p 
violates chiral symmetry. So, it is natural to separate (0|gg|0) as a factor in Eq.(17). 
The universal quark flavour independent quantity x is called magnetic susceptibility 
of quark condensate. Its numerical value had been found in [21] using a special sum 
rule: 

X = -(5.7±0.6)GeV 2 (18) 

Another example is external constant axial isovector field the interaction of which 
with light quarks is described by Lagrangian 

L' = (u'j^u - d'j fM 'j 5 d)A ll (19) 

In the presence of this field there appear induced by it condensates: 

(0|^75«|0)a = ~(0\d W \0) A = flA, (20) 

where f n = 131MeV is the constant of n — > fiu decay. The right-hand part of eq.(20) 
is obtained assuming m U) md — > 0, m 2 — > and follows directly from consideration 
of the polarization operator of axial currents U^ u (q) in the limit q — > 0, when 
nonzero contribution into H^ u (q) q ^o emerges only from one-pion intermediate state. 
The equality (20) was used to calculate the axial coupling constant in /5-decay q A 
[22]. An analogous to (20) relation holds in the case of octet axial field. Of special 
interest is the condensate induced by singlet (by flavours) constant axial field 

<0|7>> = 3/ 2 4°) (21) 

3 { S = «7m75« + hu,lhd + S7^7ss (22) 
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and Lagrangian of interaction with external field has the form 

i' = W (23) 

Constant f cannot be calculated by the method used when deriving eq.(20), since 
singlet axial current is not conserved by virtue of anomaly and the singlet pseu- 
doscalar meson rf is not Goldstone one. Constant /q is proportional to topologicaal 
susceptibility of vacuum [23] 

fo = ^Y(O), (24) 

where Nf is the number of light quarks, Nf = 3, and the topological susceptibility 
of the vacuum x(<? 2 ) is defined as 

X (q 2 ) = i J d 4 xe^(0\TQ 5 (x), Q 5 (0)|0> (25) 

Qfa) = ^g; u (x)g;„(x) (26) 

Using the QCD sum rule, one may relate /g with the part of proton spin E, carried 
by quarks in polarized ep (or /ip) scattering [23]. The value of /g was found from 
the selfconsistency condition of obtained sum rule (or from the experimental value 
of E): 

fl = (2.8 ± 0.7) • 10- 2 GeV 2 (27) 

The related to it value of the derivative at q 2 = of vacuum topological susceptibility 
x'(0), (more precisely, its nonperturbative part) is equal to: 

X'(0) = (2.3 ± 0.6) • 10~ 3 GeV 2 (28) 

The value x'(0) is of essential interest for studying properties of vacuum in QCD. 

4. Test of QCD at low energies on the basis of r-decay data. 
Determination of a s (ml) and of condensate values. 

Recently, collaborations ALEPH [24], OPAL [25] and CLEO [26] had measured 
with a good accuracy the relative probability of hadronic decays of r-lepton R T = 
B(t — + v T + hadrons) / B(t — > u T eT7 e ), the vector V and axial A spectral functions. 
Below I present the results of the theoretical analysis of these data basing on the 
operator product expansion (OPE) in QCD [27, 28]. In the perturbation theory 
series the terms up to a 3 will be taken into account, in OPE - the operators up to 
dimension 8. 

Consider the polarization operator of hadronic currents 

K = * J e^(TJ,(x)J u (0y)dx = (q,q u - q, u q 2 )Ilf (q 2 ) +q^uf(q 2 ), (29) 
where J = V, A; V^ = u^^d, A^ = u-f^ 5 d. 
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The spectral functions measured in r-decay are imaginary parts of IIj^(s) and 
uf(s), s = q 2 

vi/a^s) = 27rImU { yJ A (s + iO), a {s) = 2itImTlf\s + iO) (30) 

Functions Tly\q 2 ) and Tl^ (q 2 ) are analytical functions in the q 2 complex plane 
with a cut along the right-hand semiaxis starting from Am 2 , for Ily^g 2 ) and 9m 2 
for U^\q 2 ). Function Il^g 2 ) has kinematical pole at q 2 = 0. This is a specific 
feature of QCD following from chiral symmetry within massless u and d quarks and 
from its spontaneous violation. The kinematical pole appears due to one-pion state 
contribution into Ha(q), which has the form [27] 



P m 2 P 



Consider first the ratio of the total probability of hadronic decays of r-lepons into 
states with zero strangeness to the probability of r — > u T eV e . This ratio is given by 
the equality [29] 

B{t — > v T + hadronss=o) 
ti T ,v+A = — 



B(t — > v T ev e ) 



|2c f ds I s 

Q\V u d\ Sew / — 9 1 9 

J m A T \ m* j 



1 + 2^ I ( Vl + ai + ao)(s) - 2^U (s) 



(32) 



where |14<z| = 0.9735 ± 0.0008 is the matrix element of the Kabayashi-Maskawa 
matrix, Sew — 1.0194 ± 0.0040 is the electroweak correction [30]. Only one-pion 
state is practically contributing to the last term in (32) and it appears to be small: 

x2 2 

Ai4 0) = -24vr 2 ^i = -0.008 (33) 

Denote 

u(s) = Vl + ai + a 6 = /mfl^s) + n(l) A (s) + n^s)] = 2nImU(s) (34) 

As follows from eq.(31), H(s) has no kinematical pole, but only right-hand cut. It 
is convenient to transform the integral in eq.(32) into that over the circle of radius 
m 2 in the complex s plane [31]- [33]: 

= NVJ'J-jL | (l - (l + 2^) n W + A«<°> (35) 



Calculate first the perturbative contribution into eq.(35). To this end, use the Adler 
function D{Q 2 ): 

D{Q 2 ) = -2tt 2 ^y^I = R n an , Q 2 = - s > (36) 
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Figure 1: The applicability region of PT and OE in the complex plane s. In the 
dashed region PT + OE does not work. 



the perturbative expansion of which is known up to terms ~ a^. In MS regulariza- 
tion scheme K — K x — 1, K 2 — 1.64 [34], K 3 = 6.37 [35] for 3 flavours and for K A 
there is the estimate K 4 = 25 ± 25 [36]. The renormgroup equation yields 

da a/ . _ Q n+2 , Q 2 MQ 2 ) da 



w - ~ m = * rf'-Cm (37) 



in the MS scheme for three flavours ft = 9/4, ft = 4, ft = 10.06, ft = 47.23 
[37, 38]. Integrating over eq.(36) and using eq.(38) we get 

1 r a (Q 2 ) da 

n <« » = L D{ ^w) (38) 



Put /i 2 = m 2 and choose some (arbitrary) value a{m 2 ). With the help of eq.(37) 
one may determine then a{Q 2 ) for any Q 2 and by analytical continuation for any s 
in the complex plane. Then, calculating (38) find n(s) in the whole complex plane. 
Substitution of IT(s) into eq.(35) determines R T for the given a{m 2 ) up to power 
corrections. Thereby, knowing R T from experiment it is possible to find the corre- 
sponding to it a{m 2 ). Note, that with such an approach there is no need to expand 
the nominator in eqs.(37), (38) in the inverse powers of InQ 2 / fi 2 . Particularly, there 
is no expansion on the right-hand semiaxis in powers of the parameter n / ln(Q 2 / )i 2 ) , 
which is not small in the investigated region of Q 2 . Advantages of transformation 
of the integral over the real axis (32) in the contour integral are the following. It 
can be expected that the applicability region of the theory presented as perturbation 
theory (PT) + operator expansion (OPE) in the complex s-plane is off the shadowed 
region in Fig.l. It is evident that at positive and comparatively small s PT+OPE 
do not work. At negative s = —Q 2 in a s order a nonphysical pole appears, in higher 
orders, according with (9) it is replaced by a nonphysical cut, which starts from the 
point — Qq, determined by the formula 

, Qo f 00 da 

In ' 



- H — (39) 
Ja(u?) 13(a) 



H 2 Jafjjfl) (3(a) 

Integration over the contour allows one to obviate the dashed region in Fig.l (except 
for the vicinity of the positive semiaxis, the contribution of which, is suppressed by 
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the factor (1 — in eq.(6)), i.e. to work in the applicability region of PT+OPE. 
The OPE terms, i.e., power corrections to polarization operator, are given by the 
formula (5): 




(a s -corrections to the 1-st and 2-d terms in eq. (39) were calculated in [39] and [40] , 
respectively). Contributions of the operator with d = 2 proportional to m 2 , m 2 d 
and of the condensate 2(m u + m d )(0\qq\0) are neglected. (The latter is of an order 
of magnitude smaller than the gluonic condensate contribution). When calculating 
the d=6 term, factorization hypothesis was used. It can be readily seen that d=4 
condensates (up to small a s corrections) give no contribution into the integral over 
contour eq.(35). The contribution from the condensate (0§) may be estimated as 
|(08)| < lO~ 3 GeV 8 and appears to be negligibly small. R T y + A may be represented 
as 

R ry+A = 3^ Sew ( 1 + S' EW + 5 (0) + 5$ ] +A ) + AR® = 

= 3.475 ±0.022 (41) 

where S' em = (5/127r)o; em (m 2 ) = 0.001 is electromagnetic correction [41], = 
— (5 ± 2) • 10 -3 is the contribution of d=6 condensate (see below) and 5^ is the 
PT correction. The right-hand part presents the experimental value obtained as 
a difference between the total probability of hadronic decays R T = 3.636 ± 0.021 
[42] and the probability of decays in states with the strangeness S — — 1 R T S = 
0.161 ± 0.007 [43, 44]. For perturbative correction it follows from eq.(41) 

5 (0) = 0.206 ± 0.010 (42) 

Employing the above described method in ref.[28] the constant a s {m 2 ) was found 
from (42) 

a s (m 2 T ) = 0.355 ±0.025 (43) 

The calculation was made with the account of terms ~ a 3 ., the estimate of the effect 
of the terms ~ is accounted for in the error. May be, the error is underesti- 
mated (by 0.010-0.015), since the theoretical and experimental errors were added in 
quadratures. 

I determine now the values of condensates basing on the data [24] -[26] on spectral 
functions. It is convenient first to consider the difference Ily — 11^, which is not 
contributed by perturbative terms and there remains only the OPE contribution: 

n«M-n«( s )= £ <gL(i + CD £i) (44) 
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The gluonic condensates contribution falls out in the V — A difference and only the 
following condensates with d=4,6,8 remain 

0\~ A = 2(m u + m d ) <qq>= - /> 2 (45) 
Ol~ A = 2na s ( (u^X a d)(d^X a u) - (u l5 ^X a d)(d l5l „X a u) ) = 



= ~ —g- 1 < qq > 2 (46) 

0%~ A = 8tto s m 2 <qq> 2 , (47) 

where is determined in eq.(9). In the right-hand of (46), (47) the factorization 
hypothesis was used. Calculation of the coefficients at a s in eq.(44) gave C4 = 4/3 
[39] and c 6 = 89/48 [40]. The value of a s (m 2 T ) (43) corresponds to a s (lGeV 2 ) = 0.60. 
Thus, if we take for quark condensate at the normlization point /i 2 = lGeV 2 the 
value (6), then vacuum condensates with the account of ^-corrections appear to be 
equal (at fi 2 = lGeV 2 ): 

4 = -4.22 • 10~ 4 GeV A (48) 

6 = -3.75 • 10~ 3 GeV 6 (49) 

Og = 2.5 • 10- 3 GeV 8 (50) 

(In what follows, indeces V — A will be omitted and Ob will mean condensates with 
the account of a s corrections). 

Our aim is to compare OPE theoretical predictions with experimental data on 
V — A structure functions measured in r-decay and the values of Oq and Og found 
from experiment to compare with eqs.(49),(50). Numerical values of 6 and 8 
(49), (50) do not strongly differ. This indicates that OPE asymptotic series (44) at 
Q 2 = — s ~ lGeV 2 converge badly and, may be, even diverge and the role of higher 
dimension operators may be essential. Therefore it is necessary: either to work at 
larger Q 2 , where, however, experimental errors increase, or to improve the series 
convergence. The most plausible method is to use Borel transformation. Write for 
Ily^ — 11^ the subtractionless dispersion relation 

nS'W-nS'M^f^*^ (5i) 

Itx a Jo t — s s 

(The last term in the right-hand part is the kinematic pole contribution). Put 
s — Sq ((f) — on the upper edge of the cut) and make the Borel transformation in 
So- As a result, we get the following sum rules for the real and imaginary parts of 
(51): 



ds 2 ~ fe cos(A;0)O 



roo / g \ / g \ dg 

Jo exp (^ C0S ^J C0S (]t^ sin ^J^ 1_ai)(s) ^ = f " + ^ (_) 



2fc+2 



, , k\M 2k 



(52) 

r ( s A ■ ( s ■ j\( \( \ ds VV xfcSin^^Osfc+s 

Jo eXP li^ COS0 J Sm {M^ Sm V {Vl - ai)is) 2^ = f}^ k\ M*w 

(53) 
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0.5 1 1.5 2 0.5 1 1.5 2 



Figure 2: Eq.(52): the left-hand part is obtained basing on the experimental data, 
the shaded region corresponds to experimental errors; the right-hand part - the 
theoretical one - is represented by the dotted curve, numerical values of condensates 
are taken to be equal to the central values of eqs.(54),(55); a) = 37r/4, b) = 5n/6. 

The use of the Borel transformation along the rays in the complex plane has a 
number of advantages. The exponent index is negative at n/2 < < 3ir/2. Choose 
in the region n/2 < < ir. In this region, on one hand, the shadowed area in 
fig.l in the integrals (52), (53) is touched to a less degree, and on the other hand, 
the contribution of large s, particularly, s > m 2 , where experimental data are 
absent, is exponentially suppressed. At definite values of the contribution of some 
condensates vanishes, what may be also used. In particular, the condensate 8 does 
not contribute to (52) at = 5n/6 and to (53) at = 27r/3, while the contribution 
of Oq vanishes at = 37r/4. Finally, a well known advantage of the Borel sum 
rules is factorial suppression of higher dimension terms of OPE. Figs. 2, 3 presents 
the results of the calculations of left-hand parts of eqs.(52),(53) on the basis of the 
ALEPH [24] experimental data comparing with OPE predictions - the right-hand 
part of these equations. 

The experimental data are best described at the values [27] 

6 = -(6.8 ±2.1) • 1(T 3 GeV 6 (54) 

8 = -(7 ±4) • 1(T 3 GeV 8 (55) 

When estimating errors in (54), (55), an uncertainty of higher dimension operator 
contribution was taken into account in addition to experimental errors. (For detils 
- see [27]). 

As is seen from the figures, at these values of condensates a good agreemeent 
with experiment starts rather early - at M 2 > 0.5GeV 2 . In paper [27] the sum rules 
for the moments and the Gaussian sum rules were also considered. All of them agree 
with the values of condensates (54), (55), but the accuracy of their determination 
is worse. The values (54), (55) are by a factor of 1.5-2 larger than (49), (50). As 
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was discussed above, the accuracy of (49), (50) is of order 50%. Therefore, the most 
plausible is that the real value of condensates 6 , 8 is somewhere close to the lower 
edge of errors in (54), (55). 

Consider now the polarization operator Tl(s) defined in (34) and condensates 
entering OPE for Tl(s) (see (40)). In principle, the perturbative terms contribute 
to chirality conserving condensates. If we will follow the separation method of 
perturbative and nonperturbative contribution by introducing infrared cut-off [6, 7], 
then such a contribution would really appear due to the region of virtualities smaller 
than fi 2 . In the present paper, according to [28], an another method is exploited, 
when the /9-function is expanded only in the number of loops, (see eq.(ll) and 
the text after it) but not in 1/lnQ 2 . So, the dependence of condensates on the 
normalization point ji 2 is determined only by perturbative corrections, as is seen in 
(40). Condensates determined in such a way may be called n-loop ones (in the given 
case - 3-loop). Consider the Borel transformation of the sum Ii{s) pert + H(s) nonpert 
where U(s) per t is given by eq.(38), and H(s) nonper t - by eq.(13). Fig.4 presents the 
results of 3-loop calculation for two values of a s {m 2 ) - 0.355 and 0.330. The widths 
of the bands correspond to theoretical error taken to be equal to the last accounted 
term K 3 a 2 in the Adler function (36). (The same result for the error is obtained if 
one takes 4 loops in /5-function and puts K± = 50 ±50). The dotted line corresponds 
to the sum of contributions of gluonic condensate (11) and Oq +a condensate in (13) 
with numerical value corresponding to 0%~ A (54). The dots with errors present 
experimental data. (The contribution of the operators d = 4 and d — 6 is given 
separately in the insert). 



It is seen that the curve with a A mi 



0.330 and condensate contributions can 



be agreed with experiment, starting from M 2 = l.lGeV 2 , the agreement being im- 
proved at smaller values (0|^G 2 |0) than (11). The curve with a s {m 2 ) = 0.355 with 
the account of condensates coincides with experiment only at M 2 > 1.5GeV 2 . The 
same tendency persist for the Borel sum rules taken along the rays in the s complex 



15 



: = M 2 , GeV 2 : 

1 I .... I . I. .. I .... I .... I .... I .... I .... I . I. .. I ... . 
0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 

Figure 4: The results of the Borel transformation of V + A correlator for two values 
a s (m 2 = 0.355 and a s (m 2 ) = 0.330. The widths of the bands correspond to PT 
errors, dots with errors - experimental data. The dotted curve is the sum of the 
perturbative contribution at a s (m 2 ) = 0.330 and O4, Oq condensates. 

plane at various 0. Fig. 5 gives the sum rule for = 5n/6. From consideration of 
this and of other sum rules there follows the estimation for gluonic condensate: 

(0 | ^G%G% | 0) = 0.006 ± 0.012 GeV 4 (56) 

The best agreement of the theory with experiment in the low Q 2 region (up to 
~ 2% at M 2 > 0.8GeV 2 ) is obtained at a s {m 2 ) = 0.330 which corresponds to 
a s (m 2 z ) = 0.118. 

Let us now make some remarks on modifications of QCD in the low energy region. 

1. Analytical perturbative QCD [45], [46]. It is assumed that a s (q 2 ) is 
an analytical function of q 2 [45], or, in a more general case, it is supposed, that the 
perturbative part of the polarization operator is an analytical function of q 2 . The 
comparison of this approach with the r-decay data showed [28] that in the analytical 
QCD 

a1 na \m 2 z ) = 0.140 (57) 
what is in a strong disagreement with the world mean value a s {m 2 z ) = 0.119±0.002. 

2. Renormalon summing leading to the tachion mass A 2 in gluon propagator 
[47]. The restriction to the tachion mass 

- A 2 = 0.1 ±0.15 GeV 2 (58) 

was found from r-decay. 

3. Instantons. It was shown [28], that in the dilute instanton gas appoximation 
[48] instantons do not practically affect determination of a s (m 2 ) and the Borel sum 
rules. Their effect, however, appears to be considerable and strongly dependent on 
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1.1 1.2 1.3 1.4 1.5 



Figure 5: The result of the Borel transformation along the ray at <p = 5n/6. The 
dotted line corresponds to the central value of gluonic condensate (56), added to the 
PT curve with a s {m^) = 0.330. 



the value of the instanton radius p c in the sum rules obtained by integration over 
closed contours in the complex plane s at the radii of the contours s < 2GeV 2 . 

5. Sum rules for charmonium and gluonic condensate. 

The value of gluonic condensate had been found by Shifman, Vainstein, and 
Zakharov from the sum rules for polarized operator of vector currents of charmed 
quarks [5]. But in these calculations , the constant a s was taken comparatively 
small (a s (lGeV 2 0.3; J^qcd = WOMeV) and perturbative corrections were 
taken into account only in the first order. It is clear now, that a s (Q 2 ) in the region 
Q 2 ~ 1 4- lOGeV 2 is approximately twice as large, so that the account of higher 
order corrections became necessary. (In what follows I formulate the main results 
of [49]). 

Consider the polarization operator of charmed vector currents 

% J dx ( TJ^MO) ) = ( q»q v - g, u q 2 ) IT(g 2 ) , J, = (59) 
The dispersion representation for II (q 2 ) has the form 

R(s) = 47rImII( S + ;0) , II(g 2 ) = £ f° , (60) 

AtT z J4m 2 s{s — q A ) 

where R{po) = 1 in partonic model. In approximation of infinitely narrow widths 
of resonances R(s) can be written as sums of contributions from resonances and 
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continuum 

3 7T 

= n2 2 n ^ "ty^-ee S(s - 771 J) + 6>(s - S ) (61) 

where Q c = 2/3 is the charge of charmed quarks, so - is the continuum threshold 
(in what follows ^/sq = 4.6GeV), a(s) - is the running electromagnetic con- 
stant, a(m 2 j^) = 1/133.6 Following [5], to suppress the contribution of higher 
states and continuum we will study the polarization operator moments 

2 _ 4tt 2 ^ d \ n TJ/ f°° R(s)ds 



M »<^ s -sr \-w) n( - Q) = JL. (7f#^ (62) 

According to (61) the experimental values of moments are determind by the equality 

MniQ2) = d£ % K + 5)^ + n( So + Q2)n ^ 

It is reasonable to consider the ratios of moments M n i(<5 2 )/M n2 (Q 2 ) from which 
the uncertainty due to error in Yj/^ ee markedly falls out. Theoretical value for 
Il(g 2 ) is represented as a sum of perturbative and nonperturbative contributions. It 
is convenient to express the perturbative contribution through R(s), making use of 
(60), (62): 

R(s) = ( S ,/i 2 K(/i 2 ) (64) 



n>0 



where a(/z 2 ) = a s (fi 2 ) /ir. Nowadays, three terms of expansion in (64) are known: 
R^ [51] R^ [52], R^ [53]. They are represented as functions of quark velocity 



v = y 1 — 4m 2 / s, where m - is the pole mass of quark. Since they are cumbersome, 
I will not present them here. 

Nonperturbative contributions into polarization operator have the form (re- 
stricted by d=6 operators): 

IW ri (Q 2 ) = ^p(0 | ^G 2 | 0)[ /<«»(*) + af^(z) ] + 

+ (4mT^ /a6C(0 1 G ^ G ^ I Z = ^ ^ 

Functions f^(z), f^(z) and F(z) were calculated in [5], [54], [55], respectively. 
The use of the quark pole mass is, however, inacceptable. The matter is that in 
this case the PT corrections to moments are very large in the region of interest and 
perturbative series seems to diverge. Thus, for instance, at 

n = 10 , Q 2 = 4m 2 : ^ = 13.836 , ^ = 193.33 , = 13.791 

(66) 

(here A/W mean the coefficients at the contributions of terms ~ a to the moments, 
- are the similar coefficients for gluonic condensate contribution. In the region 
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Figure 6: The dependence of m(m) on (0\as/n)G 2 \0) obtained at n — 10, Q 2 
0.98 • 4m 2 and a s (Q 2 + m 2 ). 



of interest a ~ 0.1). At Q 2 = the situation is even worse. So, it is reasonable to 
turn to to MS mass m(/z 2 ), taken at the point fi 2 = m 2 . After turning to the MS 
mass m(m 2 ) we get instead of (66): 

M« M( 2 ) M( G 4) 

n = 10 , Q 2 = Am 2 : ^ = 0.045 , ^ = 1.136 , ^ = -1.673 

(67) 

At a ~ 0.1 and at the ratios of moments given by (67) there is a good reason to 
believe that the PT series well converges. Such a good convergence holds (at n > 5) 
only in the case of large enough Q 2 , at Q 2 = one does not succeed in finding such 
n, that perturbative corrections, a s corrections to gluonic condensates and the term 
~ (G 3 ) contribution would be simultaneously small. 

It is also necessary to choose scale - normalization point ji 2 where a s (fi 2 ) is 
taken. In (64) R{s) is a physical value and cannot depend on fi 2 . Since, however, 
we take into account in (64) only three terms, at unsuitable choice of fi 2 such /i 2 
dependence may arise due to neglected terms. At large Q 2 the natural choice is 
/i 2 = Q 2 . It can be thought that at Q 2 = the reasonable scale is /i 2 = m 2 , though 
some numerical factor is not excluded in this equality. That is why it is reasonable 
to take interpolation form 

f/ = Q 2 + m 2 , (68) 

but to check the dependence of final results on a possible factor at m 2 . Equalling 
theoretical value of some moment at fixed Q 2 (in the region where M^ and Mf> 
are small) to its experimental value one can find the dependence of m on ((a s /n)G 2 ) 
(neglecting the terms ~ (G 3 )). Such a dependence for n = 10 and Q 2 /Am 2 = 0.98 
is presented in Fig. 6. 

To fix both m and ((a s /7r)G 2 ) one should, except for moments, take their ratios. 
Fig. 7 shows the value of m obtained from the moment M w and the ratio M^jMyi 
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Figure 7: The dependence of m(m) on (0|(as/7r)G 2 |0) obtained from the moments 
(horizontal bands) and their ratios (vertical bands) at different a s . The left-hand 
figure: Q 2 = Am 2 , n = 10, M w /Mi 2 ; the right-hand figure - Q 2 = 8m 2 , n = 15, 
M 15 /M 17 . 

at Q 2 = Am 2 and from the moment M 15 and the ratio M 15 /M 17 at Q 2 = 8m 2 . The 
best values of masses of charmed quark and gluonic condensate are obtained from 
fig.7: 

m(m 2 ) = 1.275 ±0.015 GeV , (^—G 2 ^ = 0.009 ± 0.007 GeV 4 (69) 

Up to now the corrections ~ (G 3 ) were not taken into account. It appears that in 
the region of n and Q 2 used to find m and gluonic condensate they are comparatively 
small and, practically, not changing m, increase (((a s /ir)G 2 by 10 — 20% if the term 
~ (G 3 ) is estimated according to (13) at p c = 0.5 fm. 

It should be noted that improvement of the accuracy of T j/^ ee would make it 
possible to precise the value of gluonic condensate: the widths of horizontal bands 
in fig.7 are determined mainly just by this error. In particular, this, perhaps, would 
allow one to exclude the zero value of gluonic condensate, that would be extremely 
important. Unfortunately, eq.(69) does not allow one to do it for sure. Diminution 
of theoreticl errors which determine the width of vertical bands seems to be less real. 

6. Conclusion 

In this paper I compare the results of the recent precise measurements of r- 
lepton hadronic decays [24]- [26] with QCD predictions in the low energy region. 
The perturbative terms up to a 3 and the terms of the operator product expansion 
(OPE) up to d=8 were taken into account. It is shown that QCD with the account 
of OPE terms agrees with experiment up to ~ 2% at the values of the complex Borel 
parameter \M 2 \ > 0.8 — l.OGeU 2 in the left-hand semiplane of the complex plane. 
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It was found: 

1. The values of the QCD coupling constant a s (m 2 ) = 0.355 ± 0.025 from the total 
probability of r-decays and a s (m 2 ) = 0.330 from the sum rules at low energies. 
(The latter value corresponds to a s (m 2 ) = 0.118). 

2. The value of the quark condensate square (assuming factorization) 

a s (\ | 0) 2 = (2.25 ± 0.70) • 10~ 4 GeV 6 

and of quark-gluon condensate of d=8. 

3. The value of gluonic condensate: 

a) from the r-decay data: 

(0 | —G 2 | 0) = (0.006 ± 0.012) GeV 4 

7T 

b) from the sum rules for charmonium 

(0 | —G 2 | 0) = (0.009 ± 0.007) GeV 4 

It is shown that the sum rules for charmonium are in agreement with experiment 
when accounting for perturbative corrections ~ a 2 and for OPE terms proportional 
to ((a s /n)G 2 ) and to (G 3 ). 

The main conclusion is that in the range of low-energy phenomena under con- 
sideration, perturbation theory and operator expansion, i.e. the idea of vacuum 
condensates in QCD is in an excellent agreement with experiment starting from 
Q 2 ~ lGeV 2 . 

I am deeply indebted to K.N.Zyablyuk who had made the main calculations in 
papers [27, 28], the results of which I used here. 

The paper is supported by the grants CRDF RP2-2247, INTAS-2000-587 and 
RFFI 00-02-17808. 
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